Hund's coupling key role in multi-orbital correlations 
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We show how in multi-band materials, the Hund's coupling plays a crucial role in tuning the 
degree of electronic correlation. While in half-filled systems it enhances the correlations, in all other 
cases it pushes the boundary for the Mott transition at very high critical couplings. Moreover in 
weakly-hybridized non-degenerate systems the Hund's coupling plays the role of band-decoupler, 
causing a change from a collective to an individual band behavior, due to the freezing of orbital 
fluctuations. In this situation the physics is strongly dependent on individual filling and electronic 
structure of each band, and orbital-selective Mott transitions (or even a cascade of such transitions) 
are to be expected. More generally a heavy differentiation in the actual degree of correlation of 
different bands arises and the system can show both weakly and strongly correlated electrons. 

PACS numbers: 



I. INTRODUCTION 

Correlated materials are a continuous source of new 
and intriguing phenomena in condensed matter physics. 
The realizable stoichiometric structures of compounds 
partially filling 3d, 4d, 5d, 4f or 5f orbitals are virtu- 
ally endless and many of the explored ones have shown 
an extremely rich physics as a function of external pa- 
rameters such as pressure, temperature, magnetic field 
and chemical substitutions. 

This rich and varied physics is in many cases a direct 
outcome of electron-electron interactions and of the sub- 
tlety of the many-body physics that they yield. 

The "correlation strength", is the relevance of these 
many-body effects for the conduction electrons. These 
effects are indeed related to the interaction coupling 
strength in the conduction bands (but not always simply 
proportional to, as we will see in this paper), compared 
to their width (or more precisely the kinetic energy of 
the bare conduction electrons). The coupling strength is 
maximal in the mentioned orbitals in that they have a 
reduced spatial extension, which maximizes the coulomb 
interaction within the shell, yet they are filled after the 
bigger s and p shells of higher principal quantum num- 
ber, thanks to the aufbau "inversion" rule. Thus e.g. the 
4s are filled before the 3d orbitals, and determine the 
atomic positions such that the overlap of the 3d orbitals, 
and the corresponding bandwidth, is rather small. 

However these correlated orbitals are degenerate (5- 
fold for the d and 7- fold for the f) in isolated atoms 
and different behaviors arise in real solids depending on 
whether the crystal-field due to the surrounding atoms 
lifts or not this degeneracy. 

In materials where the degeneracy is totally lifted, only 
one band crosses the Fermi level and is thus partially 
filled, and a one orbital tight-binding model is often 
enough for a description of the low-energy physics. In 
this case the ratio between the inter-orbital interaction 
strength U and the bandwidth W is the control param- 



eter of the many-body effects. 

When several correlated orbitals contribute to the 
conduction bands instead, more parameters come into 
play, and a richer many-body phenomenology is to be 
expected. Besides the more complicated band struc- 
ture, the number of parameters needed to describe the 
electron-electron interaction is increased. Inter-orbital 
interactions are different from intra-orbital ones and 
Hund's rule, the tendency of electrons to distribute in 
different orbitals and in high-spin states, plays a role. 

In a multi-band tight-binding low-energy model the 
Hamiltonian reads in general: 

H = ^2 d,l„^^djrn'cr + Hint (1) 

where dl^^ creates an electron with spin a in orbital m 
on the site i and i^™ are the hopping amplitudes. The 
local interaction vertex can be written as: 

Hint = X! "-'"iT"""^ + ^ "2 ^ n.niriim' (2) 

i.7n i,m>m' 

~J ^ ^ [^2Sim- Sim' + {d^ini^d^ini^^im'^dim' \r\- h.C.) , 
i,m>m' 

where Cr — d, Tlinia — dinnjdinia^ T^iva — ^^^j'^ima-; ^im — 

lJ2aa''^lma''''T<y'dima, whcre T^a' IS the PauU matrices 
vector. U is the intra-orbital interaction, U' is the inter- 
orbital one ([/' — U — 2J holds from symmetry reasons 
on the matrix elements^ '") and J is the Hund's exchange 
coupling. J can be seen then as a measure of the depen- 
dence of the interaction strength on the electrons occu- 
pying the same or different orbitals, and on their mutual 
spin alignment. 

Among many tools for the investigation of correlated 
materials through the analysis of these low-energy mod- 
els Slave- variable mean-fields (the most known being the 
Slave Bosons, SB) ' and Dynamical Mean-Field Theory 
(DMFT) ' have collected many successes. These are non- 
perturbative local mean-fields and have been successful 
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in describing many aspects of correlated systems as band 
renormalization and mass enhancement, magnetic orders 
and the Mott transition. 

The latter is the the transition between a paramag- 
netic metal and a paramagnetic insulator because of 
the localization of conduction electrons due to interac- 
tions. It is relevant to many materials and it has also 
been invoked to explain the physics of high-temperature 
superconductors ' . The physics in the proximity of the 
Mott Transition has some well-defined features, that de- 
served the christening of a specific name: " Mottness" . 
These include strong quasiparticle mass renormalization 
and sizable spectral weight transfer from low-energy, co- 
herent, to high-energy, incoherent features, among oth- 
ers. 

The metallic phase in proximity to a Mott transition is 
a paradigm of strongly correlated physics, and its analysis 
has been performed in detail in the one-band Hubbard 
model'. 

At present there is less understanding of the Mott 
physics and in general of the effect of interactions in 
multi-band models. Early studies^' have clarified that 
the Mott transition is found at any integer filling and 
in N-band degenerate systems with bandwidth W and 
SU{2N) symmetric interaction, i.e. J = 0, the critical 
interaction strength Uc/W grows with the degeneracy N 
and is maximal at half-filling. 

The dependence of the Mott transition critical coupling 
Uc/W as a function of Hund's coupling J was studied in 
Refs.^""^-. In all these studies J is shown to reduce the 
Uc/W needed for the Mott transition, and thus to corre- 
late the electrons. However these studies are restricted to 
half-filled systems and the generality of this statement for 
the Mott transitions at all integer fillings is questionable. 

Recently, the role of Hund's coupling has been brought 
under the spotlight by the discovery of Iron-based 
superconductors ■ in which this coupling within the Fe 
orbitals is believed to be sizable. Metallic iron is known 
to have a strong Hund's coupling, leading to ferromag- 
netism, and this is expected to be true also in the Fe 
atoms of these compounds. Moreover LDA-I-DMFT cal- 
culations have shown extreme sensitivity of some aspects 
of the physics of Fe-superconductors to the Hund's cou- 
pling strength'^ In these materials 6 electrons occupy 
5 correlated bands all crossing the Fermi level, mainly of 
Fe character. The application to this case of the present 
knowledge on Hund's coupling effect is not obvious. 

One of the main results of this paper (section HI) is 
that for all integer fillings of a degenerate Hubbard model 
beside global half-filling the Hund's coupling actually in- 
creases the critical U needed for the Mott transition, 
bringing it to very large values for reasonable coupling 
strength, contrary to the previously studied half-filled 
cases. 

Many other aspects of the physics of Fe superconduc- 
tors are still subject to a debate in the scientific com- 
munity, and many competing scenarios have been put 



forward, without reaching consensus thus far. The possi- 
bility of these materials being in or in proximity to an 
Orbital Selective Mott Phase in which some localized 
electrons coexist with itinerant ones has been invoked by 
several groups ^ . The mechanism possibly leading to 
an Orbital-Selective Mott Transition- ' (OSMT) in these 
materials has been first isolated in Ref. ' and shown to 
be related to the strong Hund's coupling. 

The other main result (section IV) of this paper is 
the highlighting of Hund's coupling as the main factor 
responsible in general for the arising of OSMT in real 
non-degenerate (for the presence of a crystal-field split- 
ting or different band dispersions) weakly hybridized sys- 
tems. Overriding or accompanying bandwidth differences 
and degeneracies, Hund's coupling is seen to act as a 
band- decoupler, through the suppression of orbital fluc- 
tuations. Then the correlation character of each band is 
determined by its individual filling and structure. This 
mechanism also applies to the well known anomalous su- 
perconductor Sr2RuOi as has been recently outlined in 
Ref. ■ 



II. MODELS AND TECHNIQUES 

We will solve the model eq.(l) with Slave-spin mean- 
field (SSMF)-''' and DMFT . 

SSMF is a computationally fast slave-variable mean 
field particularly suited for studying multi-band models. 
It is extremely useful for tracing entire phase diagrams 
and yields qualitatively accurate and semi-quantitative 
results, in line with the performance of usual slave- 
bosons. But unlike the latter, which become quickly 
untractable in the multi-band case, SSMF limit the pro- 
liferation of auxiliary variables and can thus tackle i.e. 
a five band model and trace a phase diagram as a func- 
tion of several parameters. As all other slave- variable 
methods, at the mean-field level SSMF approximates 
the model with a quasi-particle effective non-interacting 
model whose parameters are self-consistently determined. 
Solving a local problem for the slave variable takes into 
account the effect of interactions and yields the quasi- 
particle effective parameters, i.e. renormalization of the 
crystal-field and of the electron mass. The latter is ex- 
pressed through the quasiparticle renormalization factor 
Zm, which ranges from 1 (non-interacting system) to 
and whose vanishing signals a Mott transition 

For supporting the SSMF results, and reporting aimed 
and quantitatively accurate results, we make use of the 
well established DMFT which is also heavier to run. 
DMFT assumes the locality of the self-energy and maps 
the (multi-orbital) lattice problem onto a (multi-orbital) 
quantum impurity problem subject to a self-consistency 
condition. The impurity problem has to be solved nu- 
merically, and we use Lanczos exact diagonalization as 
impurity solver. 

In all cases we consider non-hybridized bands, i.e. 
= t^Smm'- We then discuss how our analysis is 
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affected by a finite (small) hybridization. 

The two methods being local mean-fields, in absence of 
hybridization the k-dependence enters the problem only 
through each band dispersion. Sums over momenta can 
thus be replaced by integrals over energy weighted by 
density of states (DOS) ^(e), which we assume to be 
semicircular throughout the paper, i.e.: 

D{e)^^^l-{e/D)\ (3) 

where D = W/2 is the half-bandwidth. 

Then a model is fully specified by each band's width 
and position, the latter being set by the crystal-field split- 
ting, i.e. by the set of on-site energy levels Cm = tu- 

We solve the model in different cases: N-fold degener- 
ate model in section III, 3-band model with same band- 
with and trigonal crystal-field splitting, and with three 
different bandwidths but without crystal field in section 
IV. 

We will assume zero temperature throughout the pa- 
per. 
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Figure 1: Slave-Spin Mean-Field metal- Mott insulator critical 
coupling Uc for the N-band degenerate Hubbard model as a 
function of the ratio J /U . Only fillings such that n < N are 
shown here. The ones for n > N give identical results due 
to particle-hole symmetry. For half- filled cases (2 electrons 
in 2 bands and 3 electrons in 3 bands) the Hunds coupling 
correlates the system and reduces the critical Uc/ D. For all 
other fillings on the contrary Uc/D is strongly increased. 



III. HUNDS COUPLING AND THE MOTT 
TRANSITION IN DEGENERATE SYSTEMS 

We first study the N-fold degenerate Hubbard model 
at all possible integer fillings. All N bands have D = 1 
and no crystal-field splitting is present, i.e. = 0. 

In fig.l we show the results of the SSMF. At all inte- 
ger filling n we find, as known, that for increasing U at 
fixed J/U the renormalization factor Z„i - identical for 
all bands - decreases from unity, reaching zero beyond a 
critical Uc/D and thus signaling a Mott transition. At 
J=0 our analysis confirms that Uc is maximum for half- 
filling and that for a given filling n/N , Uc increases with 
the degeneracy N (see Ref."^). 

We then study the dependence of the Uc/D on the 
value of J. At half-filling (n = 2 electrons in = 2 
bands, or n = 3 electrons m. N = i bands in figure I) 
we find that increasing J/U reduces the value of Uc, as 
reported in previous works^""'^'''". 

But surprisingly, for any other filling Uc is rather in- 
creased by Hund's coupling. 

We see this in detail for the 2-band case in figure 2. 
Upon going from J = to a typical value J/U = 0.25, 
the critical interaction strength is reduced to a half of its 
original value at half-filling, but it increases to more than 
3 times the J = value, for the quarter-filled case. 

This result is general and does not depend on the num- 
ber of bands. In the three-band model the reduction of 
Uc at half-filling is even stronger than in the two-band 
case, while the behavior of the system with 1 electron 
follows the trend of the analogous case in the two-band 
system. Interestingly a slightly different behavior can be 
noticed in the case of 2 electrons in 3 bands. Here the on- 
set of J induces initially a reduction of Uc, which quickly 



starts increasing again and follows the trend of the other 
non-half-filled cases, with the Mott transition pushed at 
very strong couplings. 

Hints of this physics can be seen in Ref. where the 
opposite trend of the different (n — 1,2,3) Mott insu- 
lating lobes of the phase diagram as a function of J was 
reported. 

In order to understand these opposite trends we will 
follow an argument based on the analysis of the atomic 
limit"'-'"''''. 

Indeed an estimate of the Uc can be obtained follow- 
ing Mott's original image for the single-band case: at 
strong coupling the local spectral function of a Mott insu- 
lator shows features based on the atomic excitations (two 
delta- functions at a distance U from one another) , broad- 
ened by the delocalization energy which can in principle 
be calculated perturbatively in t/U. In practice these 
broadened features (called the 'Hubbard bands") have 
a width in energy of the order of the bare bandwidth 
W. Thus an estimate of the Mott gap is the distance 
of the borders of these two features, i.e. A — U ~ W. 
By the same token an estimate of the critical U for the 
Mott transition is given by the vanishing of the gap, i.e. 
Uc/W — 1, which is not far from the Ud/W ~ 1.2 ob- 
tained by full-fledged DMFT. 

The full many-body physics is more complicated , and 
in this estimates one should for example take into account 
the dependence of the delocalization energy, and thus of 
the width of the Hubbard bands, on the physical param- 
eters. In practice, at least in the Hubbard model this 
is seen to be a minor effect and only very close to the 
transition one sees sizable deviations from A — U — W. 

Transposing this argument to the multi-orbital case we 
thus want to obtain the scaling with Hund's J of the Mott 
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Figure 2: Quasiparticle renormalization factor Z for the 2- 
band Hubbard model as function of the interaction strength 
U and for different values Hund's coupling J. Upper panel: 
half-filled case. Lower panel: quarter-filled case. The effect 
of the Hund's coupling is opposite in the two cases. 



gap and of the [4 for the Mott transition in the multi- 
band Hubbard model. In Fig. 3 is depicted the atomic 
excitation spectrum for a 2-orbital atom. A change in 
the atomic spectrum will reflect directly on the distance 
between the Hubbard bands and then on the Mott gap. 
The 16 states are divided in sectors with the same number 
of particles, whose distance is set by U. Instead within 
every sector, for a given particle number, the states are 
split by J. 

The global filling of the system is set by the chemical 
potential which shift the relative position of the different 
sectors. When the system is half-filled /i is such that the 
ground state is in the 2-particle sector. For any finite J 
this is 3-fold degenerate and is the triplet of S* = 1 states. 
The gap in the atomic excitation spectrum is obtained 
summing the energy distance between the ground state 
and the lowest lying levels (i?o) in the n ~\- \ and n — \ 
sectors, i.e.: 



A" 



^0 



^0 



This quantity is readily evaluated, yielding: 
A"* = [/ -I- J. 



(4) 
(5) 
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Figure 3: Spectrum of a 2-orbital atom with Hamiltonian eq. 
2. The different sectors position in energy is as depicted for 
an half-filled system, for which the chemical potential can be 
obtained by particle-hole symmetry arguments and is ^ = 
(3C/ — 5J)/2 and the ground state lies in the n = 2 sector. 
In general the distance of this with the lowest lying state in 
the n = 3 sector is E% — Eq = —fj, + 2U — 2J whereas the 
distance with the lowest lying state in the n = 1 sector is 
Eo - Eq = fj, - U + 3J. The atomic gap is then A = U + J 
and thus is enhanced by J. 



by J. In the Mott insulator then, if we assume that the 
Hubbard bands are broadened versions of the atomic ex- 
citations the gap will also scale linearly in J, and so will 
Uc/D (with opposite sign, because a larger gap will close, 
when the Hubbard bands overlap, at a smaller U). 

Let's now replicate this argument for the quarter-filled 
case. In the corresponding atomic spectrum (Fig. 4) the 
lowest lying sector is the one with n = 1. The atomic 
gap is now: 



[/- 3J. 



(6) 



which is reduced by J. Thus the Uc for the Mott insulator 
with 1 particle (or 3) in 2 bands will get enhanced by J. 

The same argument applies in general for any number 
N of bands, to show that the Mott gap gets reduced in 
all cases but at half-filling. It is easy to calculate that in 
general: 



U-3J, 
U+{N-1)J, 



\/n^N (off half-filling), 
n = N (half-filling). 



It is independent of /i and shows that the gap is enlarged 



(7) 

The rationale behind this is the following. The energy 
of the high-spin state in each sector is lowered by the 
exchange among the electrons which is proportional to 
the number of possible pairs of aligned spare spins. The 
latter is maximum in the half-filled sector where there is 
one spare electron in every orbital. Moving away from 
half-filling it decreases quadratically, because the number 
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i ii,o>-|o,ii > of Uc with the band degeneracy The onset of J sup- 

+j presses these fluctuations rapidly and we thus observe a 

■J iti,o>+io,n> ii,i>-ii,t> quick reduction of Uc in the half-filled systems, where 

the orbital fluctuations are maximal. Moving away from 
half-filling this effect is controlled by n/N (for n < N, 
and by the particle-hole symmetric {N—n)/N for n > N) 
and it is already negligible for the system with 1 particle 
in 3 bands. The reduction of orbital correlations is in 
fact a main effect of J and we will see in section IV its 
key importance in decoupling the orbitals. 



ii,i> ii,i> + ii,t> II, i> 



Figure 4: Same as in Fig. 3 but sectors position in energy 
is depicted here for the system quarter-filled. The distance 
with the n — 2 sector is Eq — — —fi + (7 — 3 J whereas the 
distance with the n — sector is Eq — Eq — fi. The atomic 
gap is then A = U — SJ and is then reduced by J. 



of the possible pairs is Usiris — 1)/2, where Ug is the num- 
ber of spare electrons, which decreases while the number 
of either empty or fully occupied orbitals increases. 

Thus the gap in a half-filled Mott Insulator is always 
enhanced by J, being enhanced the distance between the 
ground state and the lowest state in the neighboring sec- 
tors. 

In the cases of a Mott insulator with a non half-filled 
ground state (say without loss of generality lying in the 
sector with n < N particles) this has a gain in energy due 
to J which is smaller than the one for the high-spin state 
of the neighboring sector (with n + 1 particles) closer to 
(or at) half-filling and bigger than the one with n — 1 par- 
ticles due to the mentioned quadratic dependence. The 
latter is also responsible for the fact that the subsequent 
enhancement of En-i — En with J is less than the reduc- 
tion of En+i — En- The overall gap is then reduced. 

The same obviously applies for the particle-hole sym- 
metric situation for n > N. 

The estimate of the dependence of Uc on J, following 
these arguments is: 



5Uc{n) cx 



3 J, \In^ N (off half-filling), 

-{N-1)J, n = N (half-filling). 



(8) 

In Fig. 5 we plot the SSMF phase diagram as a func- 
tion of J, and we compare the results with these asintotic 
behaviors based on the atomic limit eq. (8). It is clear 
that at strong J the atomic limit captures perfectly the 
dependence of Uc- 

However departures from these analytic results are 
found at weak J. The previous argument basis indeed, is 
that J is large enough for only the high-spin state to be 
taken into account in each sector. This assumption is not 
justified at weak J and orbital fluctuations become im- 
portant when taking into account the hopping from site 
to site. These are maximal at J = and the consequent 
gain in kinetic energy is the cause for the enhancement 
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Figure 5: (Color online) Same phase diagram as in Fig. 1, 
plotted here as a function of J. The dashed lines indicate the 
strong J behavior calculated analytically in the atomic limit, 
eq. (8) (from top to bottom SUc oc 3J, — J, — 2J). 

It is worth stressing that these conclusions are drawn 
for degenerate non-hybridized bands but are still valid 
for small hybridizations Vmm' and crystal-fields A„„„' — 
Ejn ^ £m'- Indeed both hybridization and crystal-fields 
work in general " against" J (the hybridization favors sin- 
glet states, while the crystal-field favors orbital dispro- 
portionation, both thus favor low-spin states"' ), but are 
non-singular perturbations. 

Thus if J ^ Vmm> , A„„„/ we expect these results to 
hold in general. 

More specifically one can show that their effect on the 
high-spin states in the atomic limit is null in the half-filled 
sector and lowers the energy in sectors away from half- 
filling. The detail of this energy gain depends in general 
on the specific structure of the matrices Vmm' i ^mm' and 
its analysis is beyond the scope of this paper but it can be 
easily shown for N=2 and 3 that they always reduce the 
effect of J on the gap between the atomic sectors and will 
thus reduce the effect of J on Uc- Some specific models 
with finite crystal-fields are studied in section IV. 

SSMF and atomic limit considerations are somewhat 
complementary points of view on the Mott Transition 
and they draw a quite reliable picture, but one may want 
to benchmark these results with a more quantitative and 
accurate method as DMFT. 

We plot the DMFT results for the TV = 2 degenerate 
model in Fig. C. DMFT fully confirms the SSMF scenario 
and the reliability of the SSMF method. The transition 
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boundaries are slightly shifted by DMFT as SSMF, like 
slave bosons, is known to overestimate Ur. 



in the half-filled case, thus correlating the system, while 
it lowers the self-energy in the quarter-filled case and it 
decorrelates the system. 
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Figure 6: DMFT results for 2-band degenerate Hubbard 
model for the Mott transition with n = 2 electrons (upper 
panel), and n — 1 (middle panel) which fully confirm SSMF 
results. Lower panel: phase diagram comparison between 
Slave Spin Mean-Field and Dynamical Mean-Field Theory. 
The slight overestimation of Uc is a known issue of slave- 
variable mean-fields, common i.e. to slave bosons. 

DMFT can also yield a much deeper insight into the 
physics of the systems through the analysis of dynami- 
cal quantities as i.e. the local self-energies "Smii^n) and 
spectral functions Pm{'-^) ~ —l/irlmGmi^ + ^0^^), where 
Gm is the Green Function for band m. 

We plot in Fig. 7 the self-energies at U/D = 1.75 
as a function of J for the half-filled and quarter filled 
cases. As expected J is seen to enhance the self-energy 



Figure 7: Self-energies of the metallic phase calculated in 
DMFT for U/D = 1.75 as a function of J. Upper panel: 2 
electrons in 2 bands, the self-energy grows with J. Central 
panel: 1 electrons in 2 bands, the self-energy diminishes with 
J. Lower Panel: spectral functions for the quarter-filled case 
at low energy. The features of the Hubbard bands closest to 
the broad low-energy feature (the left-most and right-most in 
the plotted energy window) move closer as J is increased, as 
predicted by the atomic limit arguments. 

In Fig. 7 we plot the spectral functions corresponding to 
the quarter-filled 2-band Hubbard model in the metallic 
phase. The Hubbard bands are very structured reflecting 
the multiple energy scales present already in the atomic 
Hamiltonian* '. In these bands, the features which are 
closer to the central peak (that is due to the coherent 
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quasi-particle excitations which disappear at the Mott 
transition) are the precursors of the lowest atomic ex- 
citations. As predicted by the atomic hmit arguments 
these features move closer to the central peak with J. 

We can thus unambiguously draw the conclusion that 
Hund's coupling works against the Mott transition in 
all cases but the half-filled one. Caveat to this state- 
ments are only the reentrant behaviour of Uc (at small 
J, where the suppression of orbital fluctuations overruns 
the shifting of the Hubbard bands) that can be found 
for integer fillings close to half, and the effect of strong 
hybridization and/or crystal field splitting in the regime 

Also, moving away from the Mott transition boundary, 
a subtle interplay of this Mott atomic-like physics with 
the Kondo screening characterizes the correlated metal- 
lic phase. In the 2-band Hubbard model these two as- 
pects work together. In the half-filled case the Kondo 
temperature is known to be reduced by J^~' '\ while for 
the quarter-filled case a modest increase of it has been 
reported^-, all in line with our results. 

For the 3-band model the situation is expected to be 
analogous both for the half-filled and the extremal n = 1 
(n = 5) cases. In the n — 2 {n — A) case instead, the 
relation between the Kondo temperature and the Uc de- 
pendence on J is probably less straightforward, in view of 
the reentrant transition boundary visible in Fig. 5 (also in 
Ref. '''' a "spin freezing transition" is reported to happen 
before the Mott transition) and this will be the subject 
of a future publication. 



IV. HUND'S COUPLING AS A 
BAND-DECOUPLER: ORBITAL-SELECTIVE 
MOTT TRANSITIONS 

Thus far we have only dealt with the fully degenerate 
Hubbard Model, where all bands are identical, and we 
have focused on the influence of the Hund's coupling on 
the behavior of the system as a whole. 

We will now allow for differences between the bands 
and study systems in presence of crystal-field splitting or 
difference between the bandwidths. We will in particu- 
larly highlight how in these systems the bands can show 
very different behaviors when J is sizable. For instance 
the bands can show a different degree of correlation, and 
the Mott transition can become selective. 

OSMT has been thoroughly studied in the 
half-filled Hubbard model with two different 
bandwidths^'*'^^'^'"'''^'^"'"' and it has been shown how 
Hund's coupling widens the parameter region where 
the OSMT can be found. However it has been shown 
in a recent publication-^ that even systems in which 
all bands have the same bandwidth can undergo an 
OSMT when J is beyond a critical value. This has been 
shown for a 3-band model populated by 4 electrons, in 
which two bands are degenerate and one is lifted up in 



energy (inset of Fig. 8). For a large range of crystal- field 
splitting the upper band can open a Mott gap becoming 
half-filled, while the other two bands remain metallic. 




U/D 



Figure 8: Phase diagram from Ref." ' for a 3-band Hubbard 
model populated by 4 electrons, with two degenerate bands 
and one lifted in energy as depicted in the inset. 

It is worth stressing the important role played by the 
Hund's coupling. As said, a finite value of J is needed 
in order to trigger the OSMT. Then further raising J/U 
widens enormously the selective phase (see Fig. 8), push- 
ing the boundary of the Mott insulating phase at very 
high U . This being a localization of 4 electrons in 3 
bands, the quick raise of the Uc necessary for this system 
to become a Mott insulator is explained by our study of 
the non half-filled Mott transition performed in section 
III. 

What is remarkable here is that while the boundary for 
the Mott transition of the whole system follows closely 
the result for the degenerate system (the result for N=3 
n=2, equivalent to n=4, of Fig. I) one band has a com- 
pletely different behavior and becomes selectively local- 
ized at a much smaller critical value of U, which gets 
further reduced by J. 

Even if this effect has been shown not very sensitive to 
a small splitting of the two degenerate bands, the main 
effect highlighted in Ref.-"^ as possibly responsible for the 
OSMT behavior is the different degeneracy of the two 
parts of the system. Albeit having the same bandwidth, 
the doubly degenerate subsystem has an enhanced kinetic 
energy that can be seen as an increased effective band- 
width (thus having a different Uc, if it were decoupled 
from the third band), tracing this OSMT back to the 
previously studied cases of systems of bands of unequal 
width. 

We will show here that the main effect causing the 
OSMT in these system of equal bands is another, namely 
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Figure 9: SSMF results for a half-filled 3-band Hubbard model 
with trigonal crystal-field splitting. One band remains half- 
filled, while the other two are shifted symmetrically, as de- 
picted in the inset. The quasiparticle weights Z\ — Z2,Zs 
are shown here, as a function of the interaction strength U/D 
for Hund's coupling J/U = 0.25 at different fixed population 
unbalance ni — n^. 



the suppression of the orbital correlations induced by J. 

For this we will study a system far from any de- 
generacy. We consider a 3-band model eq.(l) with aU 
equal bandwidths and a trigonal crystal-field splitting 
(Ai2 — — A23), at half-filling. Namely one band remains 
half-filled, one is lifted in energy, and another is symmet- 
rically lowered (inset of Fig. 9). 

Results in SSMF for a sizable Hund's couphng J/U ~ 
0.25 are shown in Fig. 9 for different fixed population 
unbalance 711 — 713. One finds that the half-filled band 
becomes insulating at a lower critical coupling than the 
other two, thus opening an Orbital-Selective Mott Phase. 
The reason for studying the model at fixed ni — 713 rather 
than at fixed A is that correlations tend to renormalize 
and strongly reduce the crystal-field splitting, so that 
keeping the population unbalance fixed the OSMT is en- 
hanced and shown to be possible. We then check the 
robustness of the OSMT at fixed A12 = 0.71? in DMFT 
and find it in a reduced but still sizable range of correla- 
tion strength (Fig. 10). 

A complete study of the phase diagram of this model is 
beyond the scope of this paper, however we have shown 
that degeneracy is not necessary for an OSMT to happen 
in a system of all equal bandwidths. 

What differentiates the bands here is only the re- 
spective individual filling. The system as a whole be- 
comes insulating at a critical U very close to the value 
= Uc{N = 3,71 = 3, J = 0.25) + A, where Uc[N = 
3,71 = 3, J = 0.25) ~ 1.151?, consistent with the atomic 
picture of the Mott insulator in which the high-spin state 
is still the ground state and the Mott Gap is reduced by 
the crystal-field splitting. 

However the half-filled band acts as "decoupled" and 



Figure 10: DMFT results for the same model as in the inset 
of Fig. 9, but fixing the crystal-field splitting A12 = — A23 = 
0.7D. Upper panel: quasiparticle residues Zm- Lower Panel: 
orbital correlations. A clear suppression of the orbital corre- 
lations between the half-filled band and the lifted ones corre- 
sponds to the OSMT. 




Figure 11: Spectral functions corresponding to the 3 phases of 
the data plotted in Fig. 10, for J = 0.25(7 and A12 = -A23 = 
0.7D. Panels from top to bottom: U/D = 1.0,1.8,2.3. The 
half-filled bands is indicated in red full line, the two shifted 
bands in blue and green dashed lines. In the top panel all 
the bands are metallic, in the bottom panel all are Mott in- 
sulating. In the middle panel only the half-filled band opens 
a Mott gap, whereas the two shifted bands are metallic, but 
the spectral density at the chemical potential violates the pin- 
ning condition imposed from the Luttinger theorem, signaling 
a non-Fermi-liquid state. 

the value of U at which it becomes insulating is smaller 
than Uc for the whole system. 

The explanation of this behavior becomes transparent 
when calculating the orbital correlations 

A^?nm' ((^7n (^77i))(^m' (^m')) 

^ {nmn„i') - {n„i) {n,n') , (9) 

in DMFT (lower panel of Fig. 10). It is clear that 
the OSMT is accompanied by an almost complete sup- 
pression of the orbital fluctuations between the half-filled 
band and the lifted ones. This implies that close to the 
transition the half-filled band is in practice decoupled 
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and it can become Mott insulating while the others re- 
main metallic. Also the effect of the crystal field will 
be minor on its Mott gap. Consequently the Uc for this 
band will be poorly affected by it. 

One may think that this suppression of the inter- 
orbital fluctuations reflects only the suppression of intra- 
orbital fluctuations typical of any Mott transition (owing 
to the fact that a Mott insulator has a very small double 
occupancy). This is not the case, here. Indeed we plot in 
Fig. 10 the ratio between these quantities 
i.e. the inter-orbital fluctuations normalized by the stan- 
dard deviation of the intra-orbital ones in each of the in- 
volved bands, which drops very fast at the OSMT. This 
shows that the inter-orbital fluctuations are suppressed 
much faster than the intra-orbital ones and highlights 
the role of this suppression driven by J in decoupling the 
bands and thus inducing the OSMT. 

In Fig. If we show the spectral functions correspond- 
ing to the three phases (metallic, orbitally-selective and 
Mott insulating) that the model goes through. The 
orbitally-selective Mott phase shows a clear gap in the 
half-fllled band only. It is also characterized by non 
Fermi-liquid behavior, here represented by a lowering of 
the density of states at the Fermi level in the metallic 
bands 

Thus one can conclude that Hund's coupling acts as a 
band decoupler, leading, to a first approximation, from a 
collective to an individual behavior for each band. Then 
it will be each band's structure and filling that determine 
the correlated properties. 

As another paradigmatic example of such a nearly de- 
coupled behavior we have studied a Hubbard model of 
3 half-filled bands with different bandwidths and Hund's 
coupling J/U = 0.25. 

We plot the results in Fig. 12. As a function of U 
the three bands undergo a "cascade" of individual Mott 
transitions, at UcS that scale very roughly with the bare 
bandwidth of each band, thus testifying an almost de- 
coupled behavior. 

The SSMF results are confirmed by DMFT as shown 
in Fig. 13, where we also plot the orbital correlations, as 
previously. The same physics is found here, with a strong 
suppression of the orbital fluctuations between localized 
and itinerant bands, leading to a decoupled behavior. 



V. CONCLUSIONS 

This work analyzes the role of Hund's coupling in tun- 
ing the correlation effects in multi-orbital materials and 
is divided in two parts. 

In the first we have shown that Hund's coupling in de- 
generate systems of interacting bands can have two oppo- 
site effects, depending on the global filling of the system: 
in half-fllled systems it strongly correlates the system and 




Figure 12: SSMF results for a 3-band Hubbard model with 
bandwidths W = 4,2,1. All bands are half-filled and J/U = 
0.25. An almost completely decoupled behavior is induced 
by J, leading to a cascade of OSMT whose critical couplings 
scale roughly with the bandwidths. 




Figure 13: DMFT results for the same model as in Fig. 12. 
Upper panel: quasiparticle residue for each band. Lower 
panel: orbital correlations, which suppressed between the lo- 
calizing band and the itinerant ones, while they remain sizable 
between the two itinerant bands. 



accordingly lowers the critical interaction strength Uc for 
the Mott transition. In all other cases the Uc is quickly 
pushed at very high values by J. This is strictly true 
when J is strong, and it can be fully understood ana- 
lyzing the atomic multiplets. When J is small instead, 
the suppression of orbital fluctuations (which is maximal 
at half-filling and decreases quickly, moving away from 
it) can dominate, enhancing correlations. This enhances 
the reduction of Uc in half-filled systems, and produces a 
reentrant behavior for fillings close to half - before cross- 
ing over to the strong J behavior -, while is negligible 
otherwise. 

The second part of this work addresses non-degenerate 
systems, where crystal field splittings or bandwidth dif- 
ferences characterize the different bands. In these sys- 
tems we have shown that the Hund's coupling acts mainly 
as a band-decoupler, through the suppression of the or- 
bital correlations. At intermediate to strong J, the bands 
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will then have different correlation effects, based on the 
individual structure or filling. When bands have simi- 
lar width and structure, the individual proximity to half- 
filling becomes then a measure of the correlation strength 
of each band. 

In systems where the two effects combine (non-half fill- 
ing of the global system and proximity to half-filling of 
one band or a subset of bands) in presence of sizable 
Hund's coupling a strong differentiation can be expected, 
because correlations in the subsystem will be enhanced 
by J whereas they can be reduced for the rest of the sys- 
tem. Strongly and weakly correlated electrons can thus 



be found in the same conduction bands. 

Materials of interest with strong Hund's coupling 
which arguably fall in this category are i.e. the Iron 
Superconductors and CaxSr2-xRuOir 
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solution disappears because the quasiparticle weight van- 
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ization factor Zm calculated in SSMF. Although Ud and 
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can be safely assumed that they both increase or they both 
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energy at low iujn. This represents an infinite quasiparticle 
lifetime at the Fermi energy which however is very sensi- 
tive to temperature. Accordingly, when using a finite low- 
energy resolution, as is done by the discrete ED impurity 
solver (9 total sites) that we used here, this strong tem- 
perature dependence shows up as a strong dependence on 
the low-energy resolution, leading strong violations of the 
pinning condition of the zero energy value of the spectral 
function of the metallic bands, as is reported in Fig. 11. 
It is worth stressing that the models studied in this work 
have zero hybridization between the bands. The conclu- 
sions of this work are indeed valid for weakly hybridized 



systems, i.e. they can be extended to nonzero hybridiza- 
tions with some care. The notion of weakly hybridized sys- 
tem is a basis-independent one. Indeed even if a change 
of basis can always diagonalize the non-interacting band 
structure, Wannier orbitals in the new basis are not local- 
ized, in general and the interaction matrix is very non-local 
then. We call a "weakly hybridized system" one that can 
be rotated in a basis where both the interaction is local 
and the hybridization between the orbitals is small. These 
two conditions cannot be satisfied for all bandstructures, 
and are on the contrary a significant constraint, which is 
material-specific . 



